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Abstract 

We extend the minimal supersymmetric standard model with bi- 
linear R-parity violation to include a pair of Higgs triplet superfields. 
The neutral components of the Higgs triplets develop small vacuum 
expectation values (VEVs) quadratic in the bilinear R-parity breaking 
parameters. In this scheme the atmospheric neutrino mass scale arises 
from bilinear R-parity breaking while for reasonable values of param- 
eters the solar neutrino mass scale is generated from the small Higgs 
triplet VEVs. We calculate neutrino masses and mixing angles in this 
model and show how the model can be tested at future colliders. The 
branching ratios of the doubly charged triplet decays are related to 
the solar neutrino angle via a simple formula. 



1 Introduction 



Ever since its proposal as a way to generate neutrino masses, either on its own 
or in the context of effective seesaw schemes pQ, many variants of the triplet 
model have been considered 1 . In fact the triplet is a common ingredient in 
the formulation of seesaw schemes with gauged [2HH1II], or ungauged B-L 
symmetry (HUE] - They are also a characteristic feature in several left-right 
symmetric models 0, required in order to reduce the extended gauge struc- 
ture down to the minimal SU(2)®U(1). For studies of the phenomenology of 
left-right symmetric models, both supersymmetric and non-sup ersymmetric 
see ref. [B] and references therein. 

Recent experiments, including the recently published first results of Kam- 
LAND have confirmed the LMA-MSW oscillation solution to the solar 
neutrino problem JU]- Together with the earlier discoveries in atmospheric 
neutrinos jTT], one can now be fairly confident that all neutrino flavours 
mix and that at least two non-zero neutrino masses exist. This impressive 
progress has brought the quest for an understanding of the smallness of neu- 
trino masses to the center of attention in particle physics. 

In the standard model neutrinos are massless. One of the best-known 
mechanisms to generate small (Majorana) neutrino masses is the seesaw 
mechanism. Elegant as it may be, the seesaw mechanism is not the only 
theoretically interesting approach to neutrino masses. More interesting from 
a phenomenological point of view are models in which the neutrino masses 
are generated at the electroweak scale. Models of this kind include, for 
example, variants of the seesaw scheme at the electro- weak scale [T^llHj. ra- 
diative models of neutrino mass [TJJE)] and supersymmetry with R-parity 
violation [T7l[THl[inil2nil2IlI221 ■ 

In this work we further explore the triplet approach to neutrino mass, 
by combining it with the idea of bilinear R-parity violating (BRpV) super- 
symmetry. The later might arise as the effective description of models with 
spontaneously broken R-parity [231121] . Alternatively, the required smallness 

^^Here the triplet majoron will not be discussed because it is ruled out by the measure- 
ment of the invisible Z decay width at LEP. 
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of the bilinear R parity violating parameters may result from some suitable 
family symmetry thus providing a common solution to the so-called 
mu-problem [2Zj and the neutrino anomalies. 

At the superpotential level lepton number violation resides only in the 
explicit BRpV terms, but after electro-weak symmetry breaking takes place, 
this induces a naturally small VEV for the neutral component of the Higgs 
triplets, suppressed by two powers of the BRpV parameters [18J, as empha- 
sized recently by Ernest Ma |2H| 2 ■ 

The result is a hybrid scheme for the neutrino masses. An attractive and 
natural possibility is that the atmospheric mass scale arises from bilinear 
R-parity breaking, while the solar mass scale is generated by the small Higgs 
triplet VEVs. The model is theoretically simple, since only tree-level physics 
is required to explain current neutrino data, and has the advantage of being 
directly testable in the next generation of colliders. 

The presence of Higgs triplets would induce lepton flavor violating decays 
of muons and taus. A number of low-energy constraints on Higgs triplet 
couplings have been derived from the non-observation of such processes ■ 
Although the model we are considering is different from the standard left- 
right models, many constraints apply equally well to our case. The most 
important bounds for our model are from the experimental upper limits on 
/it — > 3e and /i — > ej, which can be expressed as (21 



from —>■ 3e and 



Ml— 

h eu h ee < 3.2 x KT 7 - #f — — (1) 

(100 GeV) 2 V 1 

Ml- 

^ < 2 >< 10 " 6 (MteW (2) 

from fi — > cy. Here, are the Yukawa couplings of the triplet to the leptons 

and M A — is the triplet mass. Note that these bounds are on products of two 

2 This corresponds to a type II weak-scale supersymmetric seesaw scheme, where an 
effective triplet VEV is induced from the exchange of scalar bosons. This should be 
contrasted with the type-I mechanism which produces the effective neutrino mass from 
the exchange of neutralinos. Both types, however, come out proportional to the square of 
RPV parameters, as seen in eq. (tT5|l and eq. f2%| below. 
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couplings. Limits exist [2H] on individual Yukawa couplings but are weaker 
by several orders of magnitude. 

This paper is organized as follows. In the next section we will present 
the model, discussing the Higgs potential as well as the mass matrix de- 
scribing the neutrino-neutralino sector. In Sec. HU we will then turn to the 
phenomenology of the model. Neutrino masses and mixings are calculated, 
with emphasis on the solar neutrino data. Finally production and decays 
of the Higgs bosons of the model are discussed. In Sec. 13.21 we point out 
that there is a relation between various branching ratios of the decay of the 
doubly charged component of the scalar triplet and the solar angle, which 
will allow to test the validity of the model in a future accelerator. We will 
then close with a short summary. 

2 Model 

2.1 Superpotential and scalar sector 

The model to be presented below is the supersymmetric extension of the orig- 
inal triplet model of neutrino mass pQ in which the simplest form of R-parity 
violation is assumed. The model is defined by the particle content of the 
Minimal Supersymmetric extension of Standard Model (MSSM) augmented 
by a pair of Higgs triplet superfields: 
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with hypercharges Y = +2 and Y = —2 and lepton number L = —2 and 
L = +2 respectively. Note that although only one Higgs triplet superfield is 
necessary to provide the neutrinos with appropriate masses, two triplets are 
needed to avoid the triangle gauge anomaly. Apart from the new Higgs triplet 
superfields, we include R-parity violation in three generations, by adding 
bilinear lepton number violating superpotential terms [TEJ . [TT|rHHl2T)|l2 H l2~2] . 
1123 Ell- The superpotential of this model is then given by a sum of three 
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terms, 

W = Wmssm + W BRpV + W A (4) 

where 

W MS sm = ~nH u H d + h%H u Q i u c j + ^H d Q^ + h%E d t^ (5) 
is the superpotential of the MSSM, 

W B R P v = tiUHu (6) 

is the usual bilinear R-parity violating term, while 

W A = fi A A u A d + hijLitjAv (7) 

includes the terms that determine the interactions and masses of the new 
Higgs triplet superfields. The Yukawa couplings hu, hp, h B and h are 3x3 
matrices in generation space and /i, q (i — 1, . . . , 3) and M are parameters 
with units of mass. 

The scalar potential along neutral directions is a sum of two terms 

V = VsuSY + Vsoft (8) 

where 

Vsusy = \fi A A° u \ 2 + \h^u t u 3 + fi A A d \ 2 + \fiH^\ 2 + 

+ l(9 2 + 9' 2 ) {\H° d \ 2 - \K\ 2 + Ei |^| 2 + 2 |A°| 2 - 2 | A° | 2 ) 2 + 
+ Ei |2^'^A° + #°ef + \-nH° d + e L v t f 

(9) 

is the neutral part of the supersymmetric scalar potential and 

V soft = [AWuftA* - BfiH d H® + ^ + 5A/i A A° A° + c.c] + 

+M1 E, + Ml d \H»\ 2 + + 

+MiJA°| 2 + MiJA0| 2 

(10) 



4 



is the soft supersymmetry breaking scalar potential, also along neutral direc- 
tions. The expression in eq. (jlUj) contains terms which are linear in the neu- 
tral Higgs bosons H d , A°, A d and the scalar neutrinos Vi (i = 1, . . . , 3). 
The presence of bilinear RPV terms in the superpotential leads to non-zero 
sneutrino vacuum expectation values. 

The vacuum expectation values can be determined by minimizing the 
scalar potential in eq. (JHJ). These stationary conditions of the scalar potential 
are the so-called tadpole equations: 

t u = - f iBv d + ( f i 2 + M* Iu )v u + B l e l v i + V2e i h ll v3(A u ) + 

+e l €iV u - v u V 
t d = (/i 2 + M 2 Hd )v d - fi(Bv u + eVi) + v d V 
V = e i (-fiv d + B i v u + ^v j )+Mlv i + V2h»e j v u (A° u )+ 

+s/2A i ^h ij v j (A° u ) + V2fi A h^ V] (A d ) + h ij Vjh kl v kVl + (11) 

+2h^h jk v k (A° u ) 2 + v { V 
t(Ao u) = (fil + Ml u )(A° u )+v i h tJ (2h^v k (A° u ) + V2^v u ) + 

+ j= A ij h ii v . v . + B Af i A (A° d ) - 2{Al)V 

Hao) = (//| + MiJ(A°) + i M / i ^ i% + J B A/ . A (A°} + 2(A°}P 
where the non-zero vacuum expectation values are defined as 

v u = {H°), v d ={H° d ), Vi^fc) (i=,l...,3) (12) 

and 

V = \{g 2 + 9 ' 2 ){vl - vl + J>? + 2(A°) 2 - 2<A°> 2 ) (13) 

The basic idea of the model is contained in eqs. (|TT| and (|13jl and can 
be understood with the help of the following consideration. The VEVs of 
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the scalar triplets violate lepton number by two units and, therefore, are 
expected to be small. If (A°), (A°) <C v u ,Vd, 

VcV'= 1 -( g * + g'*)(vl-vl + ^) (14) 
8 i 

and one can solve eqs. (fTTj) for (A°) and (A°), 



/A°\ 1 Uij ^^[~^' J (Ml+Mi d +2X)')+J? A ^]-2^ ej ^(^+M| d +2P0 

\ A J - ^ { ^ A+ Ml u -2V>)(nl+Ml d +2V>)-(B A HA) 2 

/A°\ ^ 1 U ij fifj[-^ !J -BAMA+MA(MA+ M A u _2:D ')] _2 ^ ,: 3 1 ' u - B AMA 



(15) 



(Ml + M'l u ~ 21?') +^A d +22>') - (Bam) 2 

The Higgs triplet VEVs are quadratic in the RPV parameters, and thus 
automatically small as long as Vi -C Vd,v u . The smallness of the Vi, on 
the other hand, depends on the relative size of ej//i which we assume to be 
sufficiently smaller than one. Note also that the presence of neutral scalar 
lepton VEVs produces a mixing between scalars, higgs bosons and scalar 
neutrinos fU\ . 

2.2 Neutral Fermion Mixing 

The Lagrangian contains the following terms involving two neutral fermions 
and the neutral scalars (Higgs bosons and scalar neutrinos) 

£ D fiH®H® — eVji?° + h^UiUjA^ + h^UiUjA^ — fi A A° u A d + 



+§(M 1 A'A / + M 2 AA) + 

+^g'X'(-V i *u i - H° d *H° d + H° U *H° U + 2A°* A° - 2A°*A°)+ 



(16) 



+ _^ A 3 ( ^ V . + H 0* H _ H 0* H _ 2A 0* A + 2A 0, A } 

After the electroweak and R-parity symmetries break through non-zero vac- 
uum expectation values for the Higgs boson and sneutrinos, the Lagrangian 
contains the following mass terms involving nine neutral fermions: 

C D -^°) T M N (^°) (17) 
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where the basis is 

= K »2, -iX 3 , H° d , Hi A", AS) 

and the neutral fermion mass matrix is a 9 x 9 matrix of the form: 



(18) 



( M v 



M 



N 



IllBRpV IH//A 



m BRpV M x o 



V 



m 



m 



T 



where the symmetric matrix 

M„ = V2(A° U ) 



( h u hi2 

hl2 h-22 

hi* h 



M A 



(19) 



\ "-13 



23 



hi 3 

h>23 



(20) 



characterizes the direct contribution to the neutrino mass matrix and 
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-\g'v<i 


\g'v u 
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M 2 


\gv d 


-\gv u 






-\g'vd 


\gvd 





-n 
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\g'v u 


-\gv u 


-n 





/ 



(21) 



denotes the standard MSSM neutralino mass matrix. The matrix m B R P v is 
given as 



/ 



HlBRpV 



■\g'vi 



\gvi ei N 
■\g'v 2 \gv 2 e 2 
\g'v2, \gv z e 3 j 
and characterizes the bilinear R-parity Violation, while 



V 



(22) 



( 73^' ^ 
j-hV Vj 



(23) 



is the neutrino-Higgs triplet mixing mass matrix. On the other hand, 

/ -g>(A°) g>(A°) \ 

g(^°d) -g(K) 






m Y°A — 



V 



(24) 



/ 
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is the neutralino-Higgs triplet mixing mass matrix and 



M A 



( 



/z A 

/iA 



) 



(25) 



is the Higgs triplet mass matrix arising from the first term in eq. (|7j). 

One sees that the mixing between neutral fermions (gauginos, Higgsinos 
and neutrinos) has a rich structure, with several off-diagonal RPV entries 
which, as we show next, will induce masses (and mixings) for the neutrinos. 

3 Phenomenology 

In what follows we will work in an approximation where we neglect the 
radiatively induced neutrino masses with respect to those induced at the 
tree level due to the presence of the triplet. We have checked that there is a 
natural range of parameters in the model where the use of this approximation 
is justified. For a thorough discussion of loop-induced neutrino masses in 
BRpV see refs. jSHEOl- 

3.1 Neutrino masses and mixing angles 

One can provide an approximate analytical understanding of the tree level 
neutrino masses and mixing angles, by using the previous form of Mn. This 
form is especially convenient because the various sub-blocks in eq. (fTU|) can 
easily have vastly different orders of magnitude and one expects from eq. |T5|) 
and eq. (fTTj) 



if < jj,. More technically, in the limit where £y <C 1, where £ = m B R P v x 




■v. 



, A , m x o A < m B Rpv < M x o, M A 



(26) 




M o one can write the effective contribution of the BRpV parameters to the 
(3 x 3) neutrino mass matrix as 



M' 



eff 



eff 

Bilinear 



(27) 



8 



where the contributions from BRpV to the neutrino mass matrix is given 

by mm 

(MBiIinear) ^ " 4det(.M x o) A ^ (28) 

with 

A, = /Wi + €iV d . (29) 

With the observation eq. (J2T?|) one can further simplify the eigenvalue prob- 
lem of eq. (|27jh by means of a perturbative diagonalization. Consider first 
eq. (j2HJ). As has been discussed several times in the literature [HI], this mass 
matrix is diagonalized by only two angles, 

tanfl 13 = - Ae 1 ; (30) 

(A2+A2)2 

tanfl 23 = -^. (31) 
and leads to only one non-zero eigenvalue, given by 

where |A| 2 = X)A 2 The remaining angle 9 12 is not defined in the bilinear 
only model at tree-level. The complete diagonalization of eq. (|27|1 . M^ ff = 
R ■ M^ ff • R T , where the matrix of eigenvectors R can be expressed as a 
product of three Euler rotations, can therefore be written as, 

Motive = R . M gf Iinear ■ R T + R ■ M„ ■ R T 

= R12 • M'|?[ inear • Ri2 T + R12 • M'j, • Ri2 T 

Since Mg!f ;neaj . is already diagonal, the solar angle is defined by the entries 
of M'^ only, 

tm(2 ^> - wRi, (34) 

Note that M^, is not exactly diagonal after applying the rotation eq. (j34jl . 
However if M.' u <C Mg?[ inear , this small off-diagonals will change the angles 



9 



defined in eq. flHUj) and eq. (jSTjl only by a negligible amount. Using the exper- 
imentally measured values of tan 2 6 at m — 1 an d sin 2 26*chooz <S 1 one can 
finally find a simple formula for the solar angle (#12) which is approximately 
given by 

t fofl ^ -2y/2(h 12 - M 

tan(2d SOL ) ^ , — r — 7 777— = x (35) 

-2/iii + ft 2 2 + «33 - 2/i 2 3 

From the above discussion one expects that, for reasonable ranges of 
parameters, atmospheric neutrino physics is determined by the bilinear pa- 
rameters Aj, whereas the solar neutrino mass scale depends mostly on the 
Yukawa couplings and the triplet mass. 

Fig. Q shows an example of the solar and atmospheric neutrino mass- 
squared differences as a function of the Yukawa parameter h for a fixed value 
of the triplet mass Ma — 500 GeV. These results correspond to the following 
choice of the MSSM parameters, M 2 = 120 GeV, [i = 500GeV, tan (3 = 5, 
A = — 500GeV. In order to ensure a) negligible loop corrections due to the 
bilinear parameters and b) correct neutrino mixing angles we have chosen the 
BRpV parameters as follows: t\ = e 2 = €3 with e 2 = 10 _3 |A| and = A T = 

. Aw 

M ig 2 + M 2 g' 
4 dct ^I^O 

the numerical estimate we took the best fit Am^rpj^ = 2.5 x 10~ 3 eV 2 given 
in P2|. 

One sees that, for the hierarchical spectra produced by the model, m U2 ~ 



— ' ~ \ ' AT]VI 

10A e . The absolute value of |A| can be estimated by |A| 2 = r 2 /2 \ . In 



Amg OL scales approximately like m V2 ~ h 2 . This is expected from eq. (120 j) 
and eq. (|15j) . For values of h ~ 0(0.1), Amg OL is in the range of the LMA- 
MSW solution of the solar neutrino problem. For larger values of Ma the 
resulting m U2 gets smaller approximately like 1/Ma- From Fig. ^ one also 
sees that for large values of the Yukawas both solar and atmospheric masses 
are generated by the triplet. 

To check to which degree the simplified results for the solar angle dis- 
cussed above hold we have constructed a set of randomly chosen sample 
points and diagonalized the neutrino-neutralino mass matrix numerically. 
Points were chosen as follows. For the MSSM parameters we scan randomly 
over the following ranges: in the interval [0, 1] TeV, M 2 and \x from [0, 500] 
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0.01 0.02 0.05 0.1 0.2 0.5 1 

h 



Figure 1: Typical behavior of Am 2 ATM and Am| OL vs. the total Higgs 
triplet Yukawa coupling h, for the arbitrary choice of hu = /i23 = h/lA, 
h-22 = ^33 = hi3 = 2h/7, = 0. The triplet mass has been fixed at 
Ma = 500 GeV and the MSSM and BRpV choices are specified in the text. 

GeV, with both signs for /x, and tan/5 in [2, 15]. The resulting SUSY spec- 
tra were checked to obey existing lower limits on sparticle searches. For the 
BRpV parameters consistency with the atmospheric neutrino data requires 
|A| in the range [0.05,0.15] GeV 2 , A M ~ A T and A e < 0.3^AJ + A 2 .. 

To reduce the number of free parameters we assume Mt = //a — Ma„ = 
M& d = We then have calculated neutrino masses and mixing angles 

for several values of the triplet mass, scanning randomly over the Yukawa 
couplings with the over- restrictive constraint h = J2i<j fkj < 1- 

Fig. |21 shows the numerically calculated tan 2 ^sql versus our simple for- 
mula eq. (|35|). Clearly, for the region of interest the simple approximation 
works surprisingly well. Note, that deviations between the exact and the 
approximate results mainly occur in the region of parameter space where 
hi2 ~ h l3 . 

3.2 Implications for Accelerators 

Since in our model R-parity is violated, the lightest supersymmetric particle 
will decay. As has been shown in [HS1IS11ISH] , bilinear parameters can then 
be traced through the study of LSP decays. This feature will also remain to 
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Figure 2: Solar neutrino mixing angle vs. the ratio of triplet Yukawa cou- 
plings y ee tan 2 (^M), whe re x = z^g^gft^ and M T = 500 GeV. 

be true in the current model. We will not repeat the discussion and instead 
concentrate on the Higgs triplet in the following. 

One of the characteristic features of the triplet model of neutrino mass 
is the presence of doubly charged Higgs bosons A~~. Here we consider its 
production cross section at an e~e~ linear collider at 500 GeV center of mass 
energy [HE]- In Fig. El we present the s-channel production cross section for a 
doubly charged Higgs boson as a function of its mass. For typical expected 
luminosities of 500 inverse femtobarns per year this implies a very large 
number of events, half a million or more, for a 500 GeV mass, depending on 
the leptonic branching ratio in question. 

The next issue are the decays of such Higgs bosons. Here we come to 
the most remarkable feature of the present model, namely that the decays of 
the doubly charged Higgs bosons are a perfect tracer of the neutrino mixing 
angles. The situation here is similar to that found in the simplest bilinear 
R-parity model of neutrino mass considered in refs. |HH|l34p35j (and references 
therein). There it was found that, depending on the nature of the lightest 
supersymmetric particle, its decays patterns reflect in a simple way either 
the solar or the atmospheric mixing angles. Here we have in addition that 
the doubly charged Higgs bosons decay according to the solar mixing angle. 

In Fig.0]we give the ratio of doubly charged Higgs boson decay branching 
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100 150 200 300 500 700 1000 
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Figure 3: Cross section for the production of the Higgs triplet at center of 
mass energy equal to the triplet mass, cr{y/s = M^J, vs. the triplet mass. 
The four lines corresponds to a branching ratio (from top to bottom) of 
BR(e~e- -> A~~) = 10~\ 1(T 2 , 1(T 3 and 10" 4 . 




10" 1 10° 10 1 

tan 2 (^OL) 



Figure 4: Ratio of doubly charged Higgs boson decay branching ratios indi- 
cated by the variable y exp of eqs. (jHBj) and (|3*7)l vs. the solar neutrino mixing 
angle. The vertical band indicates currently favored values (see ref. ^U] and 
references therein). 
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ratios versus the solar neutrino mixing angle. The ratio of doubly charged 
Higgs boson decay branching ratios we consider is specified by the variable 



y exp = tan ( ] (36) 

where 



— 2y/ 2 BRn + \JIBR22 + \J2BR^ — 
with BRij denoting the measured branching ratio for the process (A~~ — ► 
lilj)- Note that the band in the plot includes an assumed 10% uncertainty 
in the measured branching ratios. The triplet mass has been fixed at Ma„ = 
500 GeV. As can be seen from the figure, there is a very strong correlation 
between the pattern of Higgs decays and the mixing angle involved in the 
solar neutrino problem. The range permitted by current solar and reactor 
neutrino data [TO] is indicated by the vertical band in Fig. |3J This correlation 
can be used as the basis for a reconstruction of neutrino angles using only 
accelerator experiments. This provides a cross-check of the determination 
provided by laboratory and underground searches for neutrino oscillations. 



4 Summary and Conclusions 

We have extended the minimal supersymmetric standard model by adding 
bilinear R-parity violation as well as a pair of Higgs triplet superfields. The 
neutral components of the Higgs triplets develop small induced vacuum ex- 
pectation values (VEVs) which depend quadratically upon the bilinear Im- 
parity breaking parameters. In this scheme, for reasonable values of pa- 
rameters, the atmospheric neutrino mass scale arises from bilinear R-parity 
breaking while the solar neutrino mass scale is generated from the small 
Higgs triplet VEVs. We have calculated the pattern of neutrino masses and 
mixing angles in this model and shown how the model can be tested at fu- 
ture colliders. The branching ratios of the doubly charged triplet decays are 
related to the solar neutrino angle via a simple formula. Similarly the atmo- 
spheric mixing can be inferred from the neutralino decay branching ratios, as 
discussed in ref. This will allow a full reconstruction of neutrino angles 
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purely from high energy accelerator experiments. The model will be tested 
in a straightforward way should a high luminosity and center-of-mass energy 
linear collider ever be built. 
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